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Abstract 


We do a precise calculation of the primordial nucleosynthesis constraint 
on the energy per length of ordinary and superconducting cosmic strings. 
For ordinary strings we provide a general formula for the constraint on 
the string tension. Using the current values for the various parameters 
that describe the evolution of loops, the constraint for ordinary strings is 
G/x < 2.2 x 10 -5 . Our constraint is weaker than previously quoted limits 
by a factor of ss 5. For superconducting loops, with currents generated 
by primordial magnetic fields, the constraint can be less stringent or more 
stringent than this limit, depending upon the strength of the magnetic 
field. We also find, in this case, that there is a negligible amount of entropy 
production if the electromagnetic radiation from strings thermalizes with 
the radiation background. 
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I. Introduction 

In the past few years there has been considerable interest in cosmic strings, both the 
ordinary and superconducting varieties. The simplest realization of these strings is 
by the breaking of a 17 ( 1 ) symmetry by a complex scalar field [l] (ordinary strings), 
and by the breaking of a 17 ( 1 ) ® U{ 1 )' symmetry by two complex scalar fields [2] 
(superconducting strings) in the early Universe. The essential parameter of the theory 
is the energy per length /i, which is related to the scale of symmetry breaking. Here we 
constrain the string tension by requiring that the inclusion of strings in our Universe 
leads to a primordial nucleosynthesis scenario that is observationally acceptable. 

Ordinary cosmic strings are cosmologically interesting because they can produce 
the density fluctuations required for galaxy formation [3], provided the energy per 
length n is G/i ss 10 -8 (a typical GUT scale) . However, loops of string radiate gravi- 
tationally at a rate that allows the energy density in gravitons to grow relative to the 
radiation background during the time that the Universe is radiation dominated. Pri- 
mordial nucleosynthesis [4] provides an upper limit to the energy density of gravitons 
produced by loops at t ss 1 sec [5,6,7], which translates to a constraint on the energy 
per length /r. Previous work [6] indicated that the constraint on the string tension 
from nucleosynthesis is a factor of a few away from ruling out galaxy formation by 
cosmic strings. 

Superconducting cosmic strings (SCS) may also explain galaxy formation, >et in 
a very different manner. Superconducting loops, with current, can radiate electro- 
magnetically as well as gravitationally. In the scenario of Ostriker, Thompson, and 
Witten [8] the currents are generated by primordial magnetic fields (pm/), and the 
electromagnetic radiation from loops of string blow bubbles in the surrounding gas, 
with galaxy formation talking place on the dense spherical shells of the gas. The suc- 
cess of their scenario depends upon the strength of the pmf . Therefore, it is of interest 
to obtain constraints on the energy per length from nucleosynthesis as a function of 
the strength of the pm/, and check for compatability with the theory of OTW. 

On the face of it, since superconducting loops also radiate electromagnetically, 
with less energy going into gravitons (compared to ordinary strings), one might think 
that the constraint on p is less stringent than that of ordinary strings. However, 
if there is a pmf that generates currents in the SCS its energy density will also 
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aifect nucleosynthesis (and the constraint on p). Another issue arises: assuming the 
electromagnetic radiation from the strings thermalizes with the radiation background, 
there will be entropy production which will dilute the energy density in gravity waves, 
loops, and the pmf relative to the background. We examine the size of this effect, 
and find that it is negligible. Also, the baryon to photon ratio rf is known at the time 
of nucleosynthesis, and at the start of matter domination. We therefore determine, 
from this knowledge, if entropy production can lead to any new constraints on the 
energy per length. 

The paper is organized as follows: in Sec II. we do a more precise analysis than has 
previously been done in obtaining a constraint on p, using primordial nucleosynthesis, 
for ordinary strings ; in Sec. III. we present a constraint on p for SCS as a fuction of 
the pmf strength ; in Sec. IV. we consider the effects of entropy production by SCS; 
in Sec. V. we summarize our work. 


II. Ordinary Strings 

We now describe a host of relations that will be essential in the forthcoming analysis. 
Loops of initial radius Lf continually form by breaking off from intersecting “infinite” 
strings when the age of the Universe is about t ss Lf/e (e 0.2, see Ref. [7]). The 
birthrate of loops of radius Lf at time t is about [9]: 

l = ™ 

(number per volume per time), where k ~ ut~ 3 ! 2 and v ss 0.01 [10] . This birthrate 
is valid from the time t , ~ t p i ane j t {Gp )~' 1 when frictional effects on the string net- 
work become negligible [11] to the time that the Universe became matter dominated. 
After formation, the number density of loops at t formed between tj and tj + dtf 
is dn(t,tf ) = KR(tf) 3 /t*jR(t) 3 dtf (assuming they have not decayed), where R is the 
FRW scale factor. For a radiation dominated Universe R <x \fi. The loops redshift 
like a nonrelativistic particle specie. 

The energy of a loop of radius L is: 


M = QpL 


( 2 . 2 ) 
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with Q = 9 1 10 ] (f3 / 2 t because the loops oscillate and are not perfectly circular). 
Loops radiate gravitationally with a power: 

Pgw = igwGp (2-3) 


where 7 gw ~ 50 — 100 [12]. 

Using the above relations, and taking into account cosmological expansion, the 
energy density at time t in gravitational waves (GW) is given by: 


GpGw(t) = 


•r(r 

J T — tm y «/ 1 j—br /1+6 


- 1.) 



m 


dr 


(2.4) 


where 9 is the unit step function, and the integrals are over the time the loops formed 
tf and the time r the radiation w r as emitted. We have also introduced the dimension- 
less parameters: 


a = *7 Gw{Gp) 2 

(2.5) 

b = iGwGpl(3e 

(2.6) 


(loops formed at t — tf will have radiated away all their energy at t = t / ( 1 + b)/b). 
These parameters are very convienent, as we do not need to specify Gp or any specific 
loop parameters. Using the standard values for the various loop parameters and 
taking Gp. ~ 10 -a yields a ~ 6 x 10~ 1J and b ~ 3 x 10 -s . At any time t, the energy 
density of loops is given by: 


Gploop 



t= * 

=f>t/l+b 


R-m) 


dt f 


(2.7) 


and the integration is over the formation times tf of the loops. A reasonable limit 
on Gp. can be obtained because of the logarithmic growth in time of pgw- Pgw ~ 
y/Upln{tjt,)!t 2 . The energy density in loops ~ y/Up/t 2 does not grow relative 
to the radiation background p ra d ~ 1 /t 2 , and for t >> t , it is small compared to pgw 
while the Universe is radiation dominated. [We shall ignore the contribution to the 
energy density in long strings, as the energy density is a factor ~ -y /Up smaller than 
that in loops.] 
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The gravitons produced by the loops are noninteracting, and will not share in the 
entropy release of other particle species. Therefore, an accurate solution to Eqns. 
(2.4) and (2.7) requires knowledge of what the rest of the Universe is composed 
of. The bulk of the energy density is in thermal equilibrium, and the equilibrium 
pressure and energy density p €q can easily be calculated by performing the ususal 
thermodynamic integrals given the masses and degeneracies of the particles (see Table 
1). Through the quark-hadron phase transition (100 MeV < T < 200 MeV) we use 
the results of Ref. [13]. The scale factor can then be calculated as a function of 
temperature: 

[. R{T)T] 3 g{T ) = constant (2.8) 


where 

sl ( 2 . 9 ) 

^ P photon 

is the effective particle degrees of freedom. In Fig. (1) we show the evolution of g{T) 
using the particles in Table (1). We see that g(T ) varies by a factor of ~ 30 between 
the GUT scale and nucleosynthesis. The evolution of the energy density in gravitons 
created from loops is paw{ t) oc g{t) 1/3 /t 2 , so the effect of the variation in g(t) is to 
decrease the expected amount of paw by roughly a half. 

To do the previous integrals we need the scale factor as a function of time. This is 
accomplished by taking steps in temperature and calculating the corresponding time 
step from the FRW equation: 

(^) J = — j ~[Peq + PGW + Plocp\ ( 2 - 10 ) 

We further take into account that neutrinos go out of equilibrium at T ~ 1 MeV, and 
do not share in the entropy release of e + , e“ annihilations. For 4 He not to be overpro- 
duced by a faster expansion rate during nucleosynthesis requires p ex tra $ 0.15p to(a j 
(at T ~ 1 MeV), and in the present case p ext ra = Pioop + Pgw- However, because 
of the time variation in ( paw + Pioop) / Ptotah we directly plug this extra contribution 
to the energy density into Wagoner’s nucleosynthesis code and require that the mass 
fraction of produced 4 He be less than 0.254. 
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We take the start of the integration time to be t . = 10 -32 sec, which corresponds 
to Gn ~ 10 -6 . To obtain an upper limit on the energy per length we must choose 
the minimum allowed values for the neutron half-life and baryon to photon ratio (at 
present):r 1 / 2 = 10.4 min., 77 = 3 x 10 _1 °. We can then determine the parameters a 
and b that give a helium mass fraction of 0.254. In Fig. (2) we plot the maximum 
allowed value for a, a max , as a function of b. The curve can be fitted, to 1%, by: 

a max = 1.54 x 10- 4 6 148 (2- 11 ) 


This relation can be used to solve for the upper limit on Gfi as a function of the loop 
parameters: 


G/xj 


= 4.39 x 10 


_0.93.-0. 02 
-8 >GW t 
^ 2 . 86 ^ 1.93 


( 2 . 12 ) 


Using the “standard” values for these parameters (7 gw = 50, (3 — 9, v — 0.01, 
e = 0.2) yields the limit: 


Gn < 2.2 x 10" 5 


(2.13) 


It must be stressed, however, that this limit is very sensitive to the parameters v and 
/?, which are somewhat uncertain. [Results from another group [14] investigating the 
evolution of cosmic strings appear to differ quantitatively with the results of Albrecht 
and Turok.] Our limit differs greatly (a factor of 5) from [6] primarily because they 
take v = 0.03. The limit is very insensitive to the start of the integration time, i.e., 
changing by « 7% for each order of magnitude we are off. Neglecting pgw in Eqn 
(2.10) results in a 23% change in the limit. 

Table (1) included only known particles, and beyond the Z boson mass scale we 
have g ^ 106.75. To investigate the effects of other possible particle species (e.g., 
supersymmetric particles) we vary g for T > 500 GeV by s x 106.75. For 3 = 
1.5,3, 10 the limit on Gfi is weaker by a factor of 1.2, 1.4, 2.1, respectively. Finally, we 
compare our limit with that obtained analytically by assuming g[T) is fixed (whence 
R(t) oc Vi) &nd taking nucleosynthesis to occur at t ss 1 sec. Using the standard 
loop parameters, our limit is a factor « 3.8 weaker than the simple estimate. Similar 
conclusions about the effect of the dilution of the gravity waves can be found in Refs. 
[6,7]. 
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In Fig. (3) we show the evolution of = pgw i Ptot and Q !oop = p loov ; p tot 

as a function of the temperature of the Universe for the choice b — 10 ~ 4 and a = 
1.84 x 10~ 10 (which results in an an extra energy density that is critically allowed). 
The occasional “bumps” in the plots represent the dilution of loops and gravity waves 
relative to the radiation background as particle species go out of equilibrium, annihi- 
late, and heat up the background. 


III. Superconducting Strings 

The procedure of the last section must now be modified to allow for the additional 
electromagnetic decay mode of the superconducting string. The electromagnetic ra- 
diation power is: 

Pem = IemJ 3 (3.1) 

where J is the current, and -jem is a numerical constant (in calculations we take 
1EM ** iGWi which is expected). [Other power laws for the EM radiation losses 
have been proposed [15], but a more thorough treatment [16] supports the form of 
Eqn. (3.1).] We use electromagnetic units here, and elsewhere, that correspond to 
e 2 = a.EMi where olem — 1/137 is the electromagnetic coupling constant. 

In the following analysis we make several restrictions. One is that the electromag- 
netic field energy will not be sufficient to stabilize the string tension, creating static 
or “floating” loops [17,18,19]. Secondly, we assume that the entropy production from 
the thermalization of the photons radiated from the SCS can be ignored, i.e., we 
assume there is no substantial dilution of gravity waves, loops, or the pmf by this 
effect. In Sec. IV we show that this assumption is appropriate. We assume that 
the pmf is a large-scale coherent field. We ignore the anisotropy introduced by the 
pmf. The effects of anisotropy on nucleosynthesis can be neglected since a measure 
of the anisotropy is the ratio of the energy density in the pmf to the total energy 
density, which is constrained by nucleosynthesis to be < 0.15. When the nucleosyn- 
thesis bound is saturated by the pmf the magnetic field strength is high enough to 
change the neutron half-life [20] and several important reaction rates. However, we 
ignore these effects which are subdominant compared to the effect of the field on the 
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expansion rate of the Universe [21]. In what follows we reduce the many degrees of 
freedom by considering only the canonical loop parameters. 

We now generalize the equations describing the energy density in loops and grav- 
itational radiation: 


Gpioop(t) 


= ~ f 


am) 


it. 


(3.2) 


Gpcw(tj = 



- 1 .) 



fi(-) 

m 


dr 


(3.3) 


where x(t,tf) = L(t,tf)/ L(tf), utf is the lifetime of a loop, and both can be deter- 
mined from the equation of motion for the loop radius. 

A superconducting loop in the presence of a pmf of strength B will have an 
induced current at formation [8]: 


Jf ss BfLf /47r/7l(u Lf) (3-4) 

where v ~ I is the expectaion value of the Higgs field. The electromagnetic radia- 
tion power for a loop formed at t/, at any time is then: 

Pem = 1 rEMpmag{tf)L*/2nln 2 (vL)L 3 (3.5) 

where the magnetic energy density p mop = B 2 /8ir. It is now useful to rewrite this as 

Pem = iGwGp. 2 f{tf) ( 3 - 6 ) 

where we have now taken iem ~ Igw and defined 

f(t f ) = p m ag{tf)L 3 f /2TrGp 3 ln 2 (vLf) (3.7) 

(£ has been approximated by Lf in the logarithm). Now p mag (tf ) oc g(tf) 1/3 /t 3 , and 
rec allin g that Lf oc tf, we see that f(tf) can vary by about two orders of magnitude 
from the time loops start forming to the onset of nucleosynthesis. Because of the time 
variation in f(tf) it becomes convienent to define a free parameter /, at the start of 
integration t 
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f(tf) = f. 


g(t,)\ U3 

9(1.)/ 


(3.8) 


which takes the place of the magnetic field. To allow direct comparison with the / of 
OTW, we note that / = f(t f ) after nucleosynthesis (after which there are no more 
particle species releasing entropy, and /(£/) becomes approximately constant). 

We assume that the drain on the primordial magnetic field by the induction of 
loop currents is negligible. A loop formed at time t will drain an energy W ^ uJj/2, 
where lj % 4nLln(vL) is the inductance of the loop. In terms of the magnetic energy 
density: 


W(t) p ma9 (f)(et) 3 //n(evt) 


(3.9) 


The evolution of the magnetic energy density is then described by: 

p mag = -4 jPma B ~ K\V(t)/t* (3.10) 

with solution: 

9™., = ftnji (3.11) 

Because ne 3 ~ 10 -3 , we cam safely treat the jrmf as noninteracting. 

The equation for the evolution of the mass of a loop is given by: 

M = 0pL = — / yGwGp. 2 [f(tf)(Lf/L) 3 + 1] (3-12) 


with the solution 

t/tf = (b+ 1 - x + y/f[tan~ 1 (x/ y/f) - fan -1 (l/ \/f)])/b (3.13) 

(the electromagnetic field energy has little effect, and has been ignored here). This 
solution is true as long as the string remains superconducting. For bosonic strings 
Jcrit ~ e y/p is the criticad current 1 , and for fermionic strings with charge carriers of 
vacuum mass M, Jcrit 55 e M. The criticad value for z, in the bosonic case is: 

*To be consistent with OTW we take this to be the critical current. However, 
the results of [19,22] indicate that J„ u = where r can be much smaller than 
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Xc Tit = Gllf(t})~lGW I <XEMlEM ~ 11.7 yjG^f(t t ) (3.14) 

(for fermionic charge carriers replace /z with A/ 3 ). The superconducting lifetime of 
loops is then: 

— = (1 + 6- + V 7[ten" 1 (W>/7) - tan-^l/ >/?)])/* (3-15) 

t f 

When x % x^t the EM field energy (~ 1/2 of the total loop energy for = e^/JI) 
will be released in the form of high energy particles. For x < x^t we are left with an 
ordinary string which will evaporate at a time t given by: 

I fcrit 

U = — = — 1 ; 

- tf b tf 

The integrals can now be performed, and Eqn (2.10) is modified to include p moB : 

r f*2ir~fGw{Gfi) 2 ln i (evt.)R(t,) A _ 157 /.(Gfi) 3 In 2 (evU)R{t.) A 

PmaB ~ 7 EM(et.) 2 R(t) A ~ R{t)Hl { ■ ’ 

which is now expressed in terms of /.. 

We restrict our attention to bosonic superconducting strings. Since the limit 
on Gfi may vary substantially in the superconducting case, we take the start of 
integration t, = tpi/(G/x) J . In Fig. (5) we show the constraint on the energy per 
length as a function of the parameter /(t/), evaluated at tf = 10 10 sec. ( roughly 
equivalent to the / of OTW). The physics is clear. For small values of / (hence 
small B fields) loops do not produce much EM radiation and the limit is the same as 
that for ordinary strings; for intermediate values of / sizeable currents are generated 
resulting in less gravitons and a weaker constraint on \i\ and larger values of / require 
smaller string tensions in order that the magnetic energy density not get too large. 
Our nucleosynthesis constraint on Gp. as a fuction of / does not restrict the (Gp., f) 
solution space of OTW that allows galaxy formation. 

unity. This can affect the OTW solution space (f,Gp) for galaxy formation, as we 
now describe. Their scenario assumes Lj > L„ it , which thereby constrains the range 
of / in their model: / < 7.3 x 10 ~ i r 2 ‘tEMhawGp. With r = 1, which they use, this 
upper limit to / is a factor ss 40 away from their solution space. Therefore, smaller 
values of r (r < 0.15) will require modification of the range of parameters that allow 
galaxy formation. 


(3.16) 
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IV. Entropy Production 

We now consider the effects of the long wavelength (~ L ) EM radiation emitted from 
the loops, which thermalize with the radiation background (during the radiation 
dominated era). Of particular interest is the time between nucleosynthesis (t as 1 
sec) and matter domination ( t zz 10 10 sec) because the entropy of the Universe is 
approximately known at these times, during which the loops might be producing an 
inconsistent amount of entropy (for a given p). Nucleosynthesis requires a baryon 
to photon ratio of 77 := (3 — 10) x 10 -10 [4]. Observations of luminous matter (i.e., 
baryons in stars) yields a lower limit of 77 ~ 3.3 x 10 -n at present. Therefore, 77 can at 
most change by a factor of % 1/30 between nucleosynthesis and matter domination. 
This information can be used to obtain a limit on Gp. 

We define = Pem to be rate of electromagnetic energy per volume produced 
by all the loops in the Universe at some temperature T. The energy released in time 
dt is then dE = T E \iR z dt, with a corresponding release of entropy dS = dE/T. In 
the radiation dominated era, 5 4pfJ 3 /3T, and the energy density p ~ 3/327T Gt 1 . 

We shall use an alternative measure of entropy production, the baryon to photon 
ratio 77 (77 oc 1 /(RT ) 3 oc 1/5). The evolution of 77, from an initial value 77; is then 
determined by: 

77/77,- = exp(— 8 ttG f T EM t 2 dt) (4.1) 

Jti 

(we ignore other forms of entropy production). 

We now calculate T E m- Looking at the energy density produced at t by loops 
formed at £/, we see that dp E M(t,tj ) = P E \t(t,tf)dni oop ,(t,tf)dt=dtfdtP E M(t,tf) 
(3 , where we have used /i(£) oc £ 1//3 . Defining y = £/£/, and letting P E m = 
7 cwGfi 3 K(y): 

GT E m = J K{y) s /ydy (4.2) 

The integral over t in Eqn. (4.1) can now be performed, 

77/77^ = ( ti / t ) s * a f KM ^ dv 


(4.3) 
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However, it would be much more useful to have x instead of y as a dependent variable, 
and from the equation of motion for a loop we solve for y in terms of i: 


y = 1 + 



dx 

1 + K(z) 


and the entropy can now be written as: 


i)/w = (fi/l)""" 1 ’'’ 


(4.4) 


(4.5) 


where the dimensionless coefficient I is obtained from 


I = 



g(f) 

1 + K{x) 



dx 

i + W) 


dx 


(4.6) 


The coefficient I was calculated for two cases: (1) tf(s) = //x 2 for x > x„ it and 
K(x) = 0 for x < x„ i t , (2) K(x) = f/x 2 for x > x„ it and K{x) = oo for i < x„ it . 
The first case corresponds to a loop with no electromagnetic radiation released when 
(and after) the critical current is reached, and the second case corresponds to all the 
loop energy going into photons at the critical current (the actual, physical situation 
should lie between these two cases). The value of 1 as a function of / is shown in Fig. 
(5) for these two cases, which give ver similar results. The coefficient I is essentially 
independent of /x, except for / >> 1 when the loops reach their critical current soon 
after they form. 

We now check our previous assumption that the pmf , loops, and gravitons will 
not have their energy density diluted by the EM radiation from loops. From Fig. (5) 
the largest possible value of G/i is ~ 10’ 4 , and from Fig. (6) the largest possible 
value of J is ~ 0.4. Using the standard loop parameters, we see that the maximum 
exponent in Eqn. (4.5) is a 3.4 x 10" 3 . For tjt = lO' 33 we have rj/rn < 0.78, and 
entropy production by the loops can be neglected. 

Requiring 17 / 17 ,' to vary by less than a factor of 1/30 over the time range t,/t = 10 ( 
yields a limit on Gfi: 


C?„ < 3.5 x 10-’-^ 

Using the standard string parameters we obtain the limit: 


(4.7) 
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Gfi <2 x 10 -2 // 2 (4.8) 

Since I < 0.4, this is far weaker than the limits obtained from nucleosynthesis. This 
is just a restatement that entropy production is negligible. 


V. Conclusion 

We have calculated the constraint on the energy per length for ordinary cosmic loops, 
and given the constraint as a function of several parameters that describe the evolution 
of loops. The bound on Gfi is clearly very sensitive to these parameters, and with 
recent studies [14,24] improving the accuracy of the parameters used here, previously 
quoted limits cannot be taken to be very accurate. 

The constraint on Gfi for superconducting strings, with small pmf, is the same 
as that for ordinary strings. Larger fields generate larger loop currents and allow 
a sizeable fraction of the loop energy to go into EM radiation, thus weakening the 
constraint. Stronger magnetic fields, with an energy density that begins to saturate 
the nucleosynthesis bound, will necessarily lead to stronger limits on fi. 

Using the bound on Gfi from nucleosynthesis we showed that there is a negligible 
amount of entropy production by SCS. Although gravitational waves emitted by loops 
are very important in constraining Gfi, the electromagnetic radiation is not. This 
situation may change, however, if the high energy, non-thermal photons resulting 
from superconducting loops making the catastrophic transition to ordinary loops is 
considered. Here, a limit is placed on Gfi by requiring that light elements are not 
overproduced by the photodisintegration of 4 He [23]. 
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Figure Captions 

1. The evolution of the particle degrees of freedom g(T). 

2. The parameter a that yields a critically allowed abundance of 4 He as a function 
of the parameter b. 

3. The evolution of flcpv and ftioop as a function of temperature for the choice 
b = 10 -4 and a = 1.84 x 10 -1 °. The abundance of 4 He in this case is critically 
allowed. 

4. The evolution of fl moJ , flew, and f hoop for the choice / = 1 and Gp = 1 x 10 -4 . 

5. The nucleosynthesis constraint for superconducting cosmic strings. We plot the 
critically allowed value of Gp as a function of /. 

6. The parameter I as a function of /. Large values of I correspond to greater 
entropy production. Case (1) corresponds to no EM radiation released at the 
critical current, case (2) corresponds to all the loop energy being released at the 
critical current. Here we have taken Gp, ss 10 - *, however, I depends strongly 
on p only for f » 1. 
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